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Abstract 

Background: We compared here the suitability and efficacy of traditional morphological approach 
and DNA barcoding to distinguish filarioid nematodes species (Nematoda, Spirurida). A reliable and 
rapid taxonomic identification of these parasites is the basis for a correct diagnosis of important 
and widespread parasitic diseases. The performance of DNA barcoding with different parameters 
was compared measuring the strength of correlation between morphological and molecular 
identification approaches. Molecular distance estimation was performed with two different 
mitochondrial markers (coxi and I2S rDNA) and different combinations of data handling were 
compared in order to provide a stronger tool for easy identification of filarioid worms. 

Results: DNA barcoding and morphology based identification of filarioid nematodes revealed high 
coherence. Despite both coxi and I2S rDNA allow to reach high-quality performances, only coxi 
revealed to be manageable. Both alignment algorithm, gaps treatment, and the criteria used to 
define the threshold value were found to affect the performance of DNA barcoding with 1 2S rDNA 
marker. Using coxi and a defined level of nucleotide divergence to delimit species boundaries, DNA 
barcoding can also be used to infer potential new species. 

Conclusion: An integrated approach allows to reach a higher discrimination power. The results 
clearly show where DNA-based and morphological identifications are consistent, and where they 
are not. The coherence between DNA-based and morphological identification for almost all the 
species examined in our work is very strong. We propose DNA barcoding as a reliable, consistent, 
and democratic tool for species discrimination in routine identification of parasitic nematodes. 


Page 1 of 12 

(page number not for citation purposes) 




Frontiers in Zoology 2009, 6:1 


http://www.frontiersinzoology.eom/content/6/1/1 


Background 

The Identification of living species is one of the major 
goals of modern biology. Species can be delimited only in 
relation to other species, it is trivial, but many discrimina- 
tore and species concepts can be used for this purpose. 
Molecular data have become widely used to aid rapid 
assessment of species diversity, and the DNA barcoding 
initiative [ 1 ] is one prominent line of research within this 
field, coordinated by the Consortium for the Barcode of 
Life (CBoL, http://barcoding.si.edu ). DNA barcoding 
involves rapid sequencing of one or a few genes from sev- 
eral representatives of a species, as well as comparisons of 
these sequences within and between species. The method 
has revealed examples of cryptic species diversity in vari- 
ous taxa [2,3]. DNA barcoding aimsatthe development of 
a universal, standardized and economical tool, but the 
fear is that to gain a sufficient accuracy the sampling 
should be massive, and, accordingly, the PCR and 
sequencing efforts expensive. Finding a balance between 
standardization, low costs and accuracy is difficult, and 
researchers have to take key decisions on the level of accu- 
racy they want to get and the costs they can reasonably 
sustain. 

A major strength of DNA barcoding is that it allows corre- 
lating any life stages of a living organism, or also a small 
part of it, to a single molecular entity (Molecular Opera- 
tional Taxonomic Unit, MOTU; sensu Floyd et al. [4] and 
does not necessarily require taxonomy skilled personnel, 
at least in the step of the molecular data generation, to be 
used in the analysis. Nevertheless, the power and pitfalls 
of the DNA barcoding approach have not yet been fully 
evaluated. In particular, the proper methods to analyse 
DNA barcoding data are stili under study (e.g. see the 
progress at CBoL working groups, http://www.bolm 
fonet.org/casestudy/index.php/browse ). 

coxi sequences are widely used for DNA barcoding of 
metazoans, but several markers have been proposed as 
putative barcodes [5,6] and different authors underline 
the importance of a DNA barcoding approach based on 
multiple markers [7,8]. Ribosomal mitochondrial genes 
are often used as alternatives to coxi marker for different 
reasons: easy to amplify, good source of synapomorphies 
in loop regions and abundance of sequences in databases 
[ 8 - 10 ]. 

Predicted problems in DNA barcoding studies include: 
poor taxonomy (e.g. single species misidentified as two or 
more species and viceversa ); insufficient sampling within a 
taxon, or insufficient sampling of taxa (see for instance 
[11] and consequent criticisms in [12-14]); polyphyletic 
or paraphyletic species [15]. 

This work focuses on an integrated approach at the Iden- 
tification of a group of nematodes, belonging to the order 


Spirurida, which includes the relevant superfamily Filari- 
oidea. Several species of filarioid nematodes are agents of 
tropical diseases both for human and other animals of 
economical value. AII the filarioids are transmitted 
through haematophagous vectore in which they span dif- 
ferent juvenile stages [16-18]. The Identification of these 
juvenile stages is a necessary condition for establishing the 
potential of transmission in endemic areas but it is diffi- 
cult, due to the small size of the juvenile stages (about 1 
mm) and paucity of charactere. Identification of juvenile 
stages is also useful to detect any possible emergent 
zoonotic filarial disease at its beginning. Laboratories typ- 
ically deals with fragments of parasitic nematodes recov- 
ered from host tissues, or with specimens representing a 
single developmental stage, and the diagnostic charactere 
are often not present in these pieces of worms. 

The identification of filarioid and related nematodes via 
DNA barcoding is an ambitious and desirable goal for 
many reasons: 1) a fast identification engine, available 
not only for taxonomists, but validated by them, is useful 
for quicker diagnoses of filariasis; 2) filarioids cause dis- 
eases of high relevance in medical and veterinary fields 
throughout the world; 3) DNA barcoding can be useful for 
those cases of difficult or impossible identification by tra- 
ditional procedures, such as co-infections with more than 
one filarioid species (e.g. Onchocerca volvulus and Loa loa ; 
see [19]); 4) parasites conferred to diagnostic laboratories 
are often of poor quality due to the difficult of sampling 
adults and undamaged organisms; 5) the model of filari- 
oid nematodes being based on a very good classical taxon- 
omy (starting from [20]) allow to avoid (as much as 
possible) problems of 'bad taxonomy' (see discussion in 
[15]); 6) DNA barcoding can offer a reliable method for 
the identification of filarioid nematodes in vectore, allow- 
ing widespread campaigns of epidemiological surveys; 7) 
nematode biodiversity is stili highly underestimated both 
at the morphological and molecular level [21], and a 
molecular approach will speed up the estimation of this 
taxonomic diversity [3]. 

Despite molecular data from representatives of filarioids 
and deposited in public databases are quite abundant for 
species of medical or veterinary relevance, very few DNA 
barcodes are available if compared with other taxa of sim- 
ilar dimensions; this is mostly caused by the difficulties of 
sampling many species of parasitic nematodes. Most of 
these DNA sequences are relative to mitochondrial genes, 
in particular 12S rDNA and coxi [9,22]. 

Here we present a double approach (morphological and 
molecular) to the taxonomic identification of filarioids 
and related nematodes on the widest (in term of species 
number) molecular collection of these parasites ever 
achieved. Morphological identification was performed by 
well known international experts, whilst molecular dis- 
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tance estimation was performed with two different rnito- 
chondrial markers ( coxi and 12S rDNA) and under 
different combinations of data handling (see below). The 
performance of DNA barcoding with different parameters 
was compared measuring the strength of correlation 
between morphological and molecular identification 
approaches. In order to provide a useful tool for easy iden- 
tification of filarioid nematodes this work aims to answer 
the following questions: 1) which is the performance of 
DNA barcoding on filarioids and related nematodes? 2) 
which is the better marker (between coxi and 12S rDNA) 
for identification these organisms at the species level? 3) 
can DNA barcoding be a useful tool for detection of puta- 
tive new species? 

Methods 

Biologica I samples, DNA extractions, PCR conditions, 

DNA sequencing and accession numbers 

Filarioids and related nematodes belong to the order 
Spirurida, a group of heteroxenic parasites with arthropod 
intermediate hosts [20,23]. In vertebrate definitive hosts, 
they are found in the digestive tract or in other different 
tissues, from the lymphatic to blood vessels and heart 
chambers, from abdominal and thoracic cavities to skin 
and subcutaneous tissues. Well preserved biological sam- 
ples are not easily obtained for these parasites, since dis- 
section of vertebrate hosts is generally required for 
collection. We emphasize that most of the specimens for 
which we generated DNA sequences derive from wild nat- 
urally infected hosts, and most of the samples have been 
recovered at necropsy. In spite of these difficulties, we 
have included in this study the most important filarioid 
parasites of humans and other animals, including 
Onchocerca volvulus, agent of human river blindness, 
Wuchereria bancrofti and Brugia malayi, agents of human 
tropical elephantiasis, Loa loa, agent of human ocular 
filariasis, Dirofilaria immitis, agent of heartworm disease of 
dogs and cats, plus a collection of specimens recovered 
from the tissues of wild animals such as bats, ungulates, 
monkeys, tropical toads, reptiles and birds, collected all 
around the world (see Table 1 for a summary of the spe- 
cies considered in this study; for further details on these 
organisms see additional file 1: 'Investigated specimens'). 

All the biological material analysed have been stored fol- 
lowing the procedures specified in the Biorepositories ini- 
tiative http://www.biorepositories.org and belong to the 
collection identified as ':zpl' of MIB institution (which 
represents a confirmed record at Biorepositories initia- 
tive). Details on parasite species included in this work are 
given in additional file 1: 'Investigated specimens' 

DNA extraction procedures, PCR conditions and sequenc- 
ing of amplified DNA fragments were performed follow- 


ing Standard procedures (details are found in the 
additional file 2: 'Experimental conditions'). Primers used 
for amplification are: coxi: coIintF and coIintR [22]; 12S 
rDNA: 12SF and 12SR [9]. The sequences generated have 
been deposited in the EMBL Data Library according to the 
EBI Barcoding Procedure (see details available at http:// 
www.ebi.ac.uk/embl/Documentation/ 

Barcode submission.html ) under the following acces- 
sions: coxi: [GenBank: AM749226 - AM749298 . Gen- 

Bank: AM886173 ] ; 12S rDNA: [GenBank: AM77976.9 - 
AM779855 ]. The detailed list of accession numbers is 
found in the additional file 1: 'Investigated specimens'. 

The datasets have been deposited in EMBL data alignment 
under the following accessions: coxi alignment [Gen- 
Bank ALIGN 001178 ]: 12S rDNA alignment: [Gen- 

Bank: ALIGN 001179 ] . 

The morphological identification procedure 

For species identification, a morphological anatomical 
analysis is performed with worms cleared in lactophenol 
and using an optical microscope equipped with a camera 
lucida. The characters studied have been validated since 
years [24] and are the basis of the key of identification 
[20]. They include the measurements, the number and 
disposition of the sensory papillae on head and male tail, 
the different parts of the digestive tract and of the genital 
apparatus. A series of other characters have been intro- 
duced for precise identification; these are thought to be 
important during mating and able to discriminate close 
species: the cuticular ornamentation of male posterior 
region, or area rugosa, which acts as anti-slit system; the 
spicule distal extremities; the muscular-hypodermal anat- 
omy. In filarioids, the first stage larva or microfilaria is a 
good discriminative character and is particularly studied: 
specimens are fixed in extension and measured; the 
cephalic hook and other cuticular head ornamentation 
are analysed as well as caudal extremity. For the correct 
observation of many characters manipulations are neces- 
sary: dissection of spicules and ovijector, head cut and ori- 
entation in front view, etc. (see [25-27]). 

Deflnition of molecular datasets 

The DNA sequences used in this study were obtained by 
direct sequencing of PCR products or collected from Gen- 
Bank; only sequences meeting a priory defined criteria of 
length, position, similarity and taxonomy were analysed. 
Each DNA sequence analysed belong to one of the four 
following groups: (1) sequences originated from organ- 
isms morphologically identified by international experts 
of our group; (2) sequences collected from GenBank and 
morphologically identified by international experts not 
affiliated to our group; (3) sequences originated from 
organisms collected by our group and morphologically 
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Table I: List of the species included in this study. List of nematodes species included in this study and their relevance in human (H), 
veterinary (V) or zoonotic (Z) parasitic diseases. Species used as models in researches are also indicated (model). 


Species 

Relevance 

Acanthocheilonema reconditum (Grassi, 1 890) 

V 

Acanthocheilonema viteae (Krepkogorskaya, 1933) 

model 

Brugia malayi (Brug, 1 927) 

H 

Brugia pahangi (Buckley & Edeson, 1 956) 
Cercopithifilaria bulboidea Uni & Bain, 2001 
Cercopithifilaria crassa Uni, Bain & Takaoka, 2002 
Cercopithifilaria japonica (Uni, 1983) 

Cercopithifilaria longa Uni, Bain & Takaoka, 2002 
Cercopithifilaria minuta Uni & Bain 200 1 
Cercopithifilaria multicauda Uni & Bain, 2001 
Cercopithifilaria roussilhoni Bain, Petit & Chabaud, 1986 
Cercopithiplaria shohoi Uni, Suzuki & Katsumi, 1998 
Cercopithifilaria tumidicervicata Uni & Bain, 2001 
Dipetalonema gracile (Rudolphi, 1 809) 

V, model 

Diroflaria (Dirofilaria) immitis (Leidy, 1 856) 

V, Z, model 

Dirofilaria (Nochtiella) repens Railliet & Henry, 1911 
Filaria martis Gmelin, 1 790 
Foleyella furcata (Linstow, 1 899) 

Litomosa westi (Gardner & Smith, 1986) 
Litomosoides brasiliensis Lins de Almeida, 1936 
Litomosoides galizai Bain, Petit, Diagne, 1989 
Litomosoides hamletti Sandground, 1934 
Litomosoides scotti Forrester & Kinsella, 1973 

V, Z 

Litomosoides sigmodontis Chandler, 1 93 1 
Litomosoides yutajensis Guerrero, Martin & Bain, 2003 

model 

Loa loa (Cobbold, 1864) 

Loxodontoflaria caprini Uni & Bain, 2006 

Mansonella (Cutifilaria) perforata Uni, Bain & Takaoka, 2004 

Mansonella (Tetrapetalonema) atelensis amazonae n. subsp. Bain & Guerrero, 2008 

Ochoterenella sp. sensu Casiraghi et al„ 2004 

H, model 

Onchocerca dewittei japonica Uni, Bain & Takaoka, 200 1 
Onchocerca eberhardi Uni & Bain, 2007 

Z 

Onchocerca gibsoni (Cleland & Johnston, 1910) 

V, model 

Onchocerca lupi Rodonaja, 1 967 

V 

Onchocerca ochengi Bwangamoi, 1 969 
Onchocerca skrjabini Ruklyadev, 1 964 
Onchocerca suzukii Yagi, Bain & Shoho, 1994 

V, model 

Onchocerca volvulus (Leuckart, 1 893) 
Piratuba scaffi Bain, 1 974 

H 

Setaria digitata (Linstow, 1 906) 

V, Z 

Setaria equina (Abildgaard, 1789) 

V 

Setaria labiatopapillosa (Alessandrini, 1848) 
Setaria tundra Issaitshikoff & Rajewskaya, 1928 

V 

Spirocerca lupi (Rudolphi, 1 809) 

V 

Thelazia callipaeda Railliet & Henry, 1910 

V 

Thelazia gulosa (Railliet & Henry, 1910) 

V 

Thelazia lacrymalis (Gurlt, 1831) 

V 

Wuchereria bancrofti (Cobbold, 1 877) 

H 


undetermined; (4) sequences collected from GenBank 
and whose Identification process is not certainly based on 
morphology. 

DNA sequences were partitioned in three types of datasets 
(called here A, B and C) based on the analyses to be per- 
formed. In order to carry out DNA barcoding study with 


the Standard marker coxi two datasets called A and B have 
been produced. 

Dataset A encompasses only sequences derived from spec- 
imens for which morphological identification was sure 
(cases 1 and 2) and was used to test the coherence 
between morphological and molecular approaches fol- 
lowing a 'classica!' DNA barcoding analysis: generation of 
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a Kimura 2-parameters (K2P; [28]) distance graph and 
cumulative error plots. coxi dataset A includes 151 
sequences 627 bp long representing 46 morpho-species 
(with an average of 3.3 specimens per species; Standard 
deviation 3.4; range: 1-20). 

Dataset B encompasses ali coxi sequences available (cases 
1, 2, 3 and 4) and was used for Standard DNA barcoding 
analyses with the most comprehensive dataset (we under- 
lined that this dataset contains also sequences derived 
from morphologically undetermined organisms. coxi 
dataset B includes 168 sequences 630 bp long (gaps are 
taken into account). 

Finally, two datasets identified as as type C (one relative 
to coxi and one relative to 12S rDNA) encompass 
sequences derived from organisms belonging to cases 1 
and 2 and for which both genes were available. These two 
types of C datasets were used to compare the molecular 
Identification performance of different markers, and of 
different data handling. The two datasets C include 86 
sequences ( coxi are 627 bp long; 12S rDNA are 643 bp 
long including gaps) representing 44 morpho-species 
(with an average of 2.0 specimens per species; Standard 
deviation 1.5; range: 1-6). 

DNA barcoding analyses 

In order to evaluate the performance of the DNA barcod- 
ing approach performed on filarioid nematodes, the 
degree of correlation between the species identification 
based on morphology and on molecular divergences was 
measured. This test was developed for the mitochondrial 
gene coxi on the widest molecular dataset of filarioid nem- 
atodes identified by morphological experts (dataset A). 

Intraspecific, interspecific, overall mean K2P distances 
[28] and relative Standard errors were calculated with 
MEGA4.1 [29] - options = Kimura 2-parameters, pairwise 
deletion. 

Typical DNA barcoding analyses are based on the compar- 
ison between intraspecific and interspecific distribution of 
nucleotide divergence that allow the inference of a molec- 
ular threshold to help taxonomic decision. Based on this 
approach two kinds of error can occur. Type I errors (false 
positive) occur when co-specific specimens show a genetic 
distance greater than threshold value. In contrast, type II 
errors (false negative) occur when genetic distance minor 
to the threshold value is found between different species. 
Cumulative error plots show the error rates generated by 
both type I and type II errors based on different values of 
threshold [ 14] . In this context, the threshold value relative 
to the minor rate of cumulative error is called optimum 


threshold (OT). When not a single value, but a range of 
threshold values is relative to the same minimum cumu- 
lative error, the formal OT is calculated as the average 
value of the range. Differently, a Standard threshold (ST) 
value was calculated as 10 times the mean intraspecific 
variability according to Hebert et al. [11]. Cumulative 
error rates relative to ST and OT were also compared. 

DNA taxonom y analyses 

According to Lefebure et al. [8] we will refer to the terms 
DNA barcoding and DNA taxonomy respectively for: 1) 
identification of organisms based on DNA sequence vari- 
ability and assignment to a certain species previously 
described; 2) prediction and classification of new taxa 
using DNA. 

On the bases of the results obtained with dataset A, the OT 
generated has been used to perform DNA barcoding and 
DNA taxonomy approaches on dataset B. 

The resulting K2P distance matrix has been used: 1) to 
infer MOTUs delimited by OT; 2) to analyse the MOTU 
composition testing the congruence with previously 
described species (DNA barcoding); 3) to perform predic- 
tion and classification of potentially new taxa (DNA tax- 
onomy). 

A phenetic tree was also generated for type B dataset of 
coxi marker using MEGA 4.0 [29] - options = tree infer- 
ence method: neighbor-joining; phylogeny test and 
options: bootstrap (100 replicates); gaps/missing data: 
pairwise deletion; codon positions: lst+2nd+3rd+non- 
coding; substitution model: K2P; substitutions to include: 
transitions + transversions; pattern among lineages: same 
(homogeneous); rates among sites: uniform rates. 
MOTUs previously identified from K2P distance matrix 
have been showed on the tree with squared brackets. 

Differentia I performance of DNA barcoding 

Different works show the importance of a proper data 
management (from the choice of alignment Software to 
the gap treatment) in the context of DNA barcoding anal- 
yses (see for instance [30]). The relevance of a DNA bar- 
coding approach based on multiple marker is also 
underlined by different authors [7,8]. 

In this work we compared the performance of DNA bar- 
coding performed with different combination of data 
handling and with different DNA barcodes. The perform- 
ance comparison was based on the measure of the 
strength of correlation between morphological and 
molecular approaches (cumulative error rates were com- 
pared). 
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In this connection homologous DNA sequences of coxi 
and 12S rDNA type C datasets were aligned with two dif- 
ferent multiple alignment Software: MUSCLE [31] - 
default options - and ClustalX [32] - default options. 
Alignments were hand corrected with BioEdit [33] in 
order to discard the terminal misalignments and were 
pruned to 627 bp for coxi and 643 bp for 12S rDNA (gaps 
included). K2P distances [28] were calculated with two 
different applications, and the gaps were treated in two 
different ways: MEGA [29] - options = Kimura 2-parame- 
ters, both pairwise deletion and complete deletion were 
set in separate runs - and TREECON [34] - options = 
Kimura 2-parameters, both 'not take into account' and 
'take into account' were set in separate runs. K2P distance 
graph and cumulative error plots were produced for the 
sixteen combination of software/parameters tested on two 
type C datasets. A schematic representation of the differ- 
ent approaches used is illustrated in Table 2. Intraspecific, 
interspecific, overall mean K2P distances [28] and relative 
Standard errors were calculated for coxi and 12S rDNA 
datasets (after alignment with MUSCLE) with MEGA 4.1 
[29] - options = Kimura 2-parameters, pairwise deletion. 

Results 

Morphologica I 'Identification 

A total of 89 specimens collected from 21 localities have 
been analysed by morphological experts. 76 specimens 
have been identified as 28 morpho-species, belonging to 
12 genera, while 1 1 specimens (forming 5 distinguishable 
morphological groups) have not been assigned, at this 
level of the work, to any described species. 

Analyses on the molecular datasets 

The datasets generated in this work comprise a total of 
254 gene sequences, 141 of which were produced in this 


study (for details see additional file 1: 'Investigated speci- 
mens'). For a total of 11 morphologically identified spe- 
cies, the DNA gene sequences here reported represent the 
first entries in GenBank. 

DNA barcoding: coherence between molecular and 
morphological identifications 

The multiple alignment of coxi gene sequences forming 
dataset A presents no insertion/deletion (indels). coxi 
mean nucleotide distance within species is 0.5% (Stand- 
ard error: 0.6%; range: 0 - 2.4%); coxi mean nucleotide 
distance between species is 16.2% (Standard error: 3.7%; 
range: 0 - 27.8%); coxi overall mean diversity is 16.0% 
(Standard error: 1.0%). 

Figure 1 shows the frequency distribution of intraspecific 
and interspecific genetic divergences in coxi dataset A. An 
overlap between the two distributions is observable at val- 
ues minor to 2%. Since some interspecific divergences are 
as low as 0% it is not possible to set any threshold value 
that allow to exclude false negatives (type II errors). 

The minimum cumulative error is 0.62% (see Figure 1) at 
a threshold level of 4.8% (OT). ST (10 times intraspecific 
mean divergence) assumes the value of 5.0% and gener- 
ates the same cumulative error (0.62%). As shown in K2P 
distance graph (Figure 2), using OT or ST, no overlap of 
intraspecific and interspecific nucleotide divergence 
occurs at distance values greater than threshold values 
(hence no false positive occur, type I errors). In contrast, 
as stated before, at distance values lower than OT or ST, a 
degree of overlap is observable (false negatives, type II 
errors). This percentage of false negatives (that represents 
the overall amount of cumulative error) are generated by 
two couples of congeneric species: 1) O. volvulus and O. 


Table 2: Minimum cumulative errors (MCE). Minimum cumulative errors relative to Standard threshold (MCE ST ) and optimum 
threshold (MCE OT ) for different markers and different data handiing. 


marker 

alignment 

distance calculation 

gap treatment* 

ST 

MCE st (%) 

OT 

MCE ot (%) 

coxi 

MUSCLE 

TREECON 

TA 

4.6 

0.3 

4.5 

0.3 

coxi 

MUSCLE 

TREECON 

NTA 

4.6 

0.3 

4.5 

0.3 

coxi 

MUSCLE 

MEGA 

PD 

4.6 

0.3 

3.9 

0.3 

coxi 

MUSCLE 

MEGA 

CD 

4.2 

0.3 

4.5 

0.3 

coxi 

Clustal 

TREECON 

TA 

4.6 

0.3 

4.5 

0.3 

coxi 

Clustal 

TREECON 

NTA 

4.6 

0.3 

4.5 

0.3 

coxi 

Clustal 

MEGA 

PD 

4.6 

0.3 

3.9 

0.3 

coxi 

Clustal 

MEGA 

CD 

4.2 

0.3 

4.5 

0.3 

I2S rDNA 

MUSCLE 

TREECON 

TA 

34.7 

99.9 

9.0 

0.4 

I2S rDNA 

MUSCLE 

TREECON 

NTA 

22.0 

88.7 

6.4 

0.3 

I2S rDNA 

MUSCLE 

MEGA 

PD 

22.0 

89.9 

6.7 

0.3 

I2S rDNA 

MUSCLE 

MEGA 

CD 
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TREECON 
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22.0 
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NTA 
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0.4 

I2S rDNA 
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MEGA 

PD 

14.3 

14.8 

5.8 

0.4 

I2S rDNA 
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MEGA 

CD 

8.1 

16.2 

4.4 

l.l 


*NTA is for not taken into account; TA is for taken into account; PD is for pairwise deletion, CD is for complete deletion 
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Figure I 

Cumulative error plot. Type I (yellow) and type II (red) errors obtained with different thresholds for coxi sequences of 46 
spirurida species. 


ochengi (mean interspecific distance is 1.9%); 2) C. bulboi- 
dea and C. longa (mean interspecific distance is 0.2%). If 
O. volvulus and C. bulboidea are discarded from dataset, no 
overlap between intraspecific and interspecific distribu- 
tions are observable, and the OT allows to reach 0% of 
cumulative error. In summary, Identification based on 
molecular divergence threshold for coxi is coherent with 
morphological approach for 44 species out of 46. 

DNA taxonomy: a direct applicatiori of DNA barcoding 

Sequences in dataset B (that encompasses all of the coxi 
sequences of filarioid nematodes provided by GenBank) 
were used to generate a new K2P distance matrix. This has 
been used to identify the MOTUs whose boundaries are 
delimited by the OT value (4.8%) calculated at previous 
point. This approach was used to reach two different pur- 
poses: DNA barcoding (MOTUs composition were 
checked to correspond to previously identified species) 
and DNA taxonomy (which allow to identify potentially 
new species). 51 MOTUs were identified: 46 belong to 
species previously described, and five belong to not mor- 
phologically identified organisms (potentially non 
described species) of filarioid nematodes. 

Molecular cryptic species (complexes composed by differ- 
ent morpho-species grouped into single MOTUs) are stili 
present: all of the sequences of C. bulboidea and C. longa 


form a single MOTU and all of the sequences of O. volvulus 
and O. ochengi form also a single MOTU. This resuit is 
concordant with the results obtained from dataset A. 

The five unidentified MOTUs encompass specimens col- 
lected from avian and mammals hosts. Three out of these 
contain specimens collected from African mammals ( Oiyx 
gazella, Redunca fulvorufula and Equus zebra hartmannae ) 
and the remaining two contain specimens collected from 
three bird species ( Sitta europea, Paradoxornis webbianus 
and Sturnus vulgaris). Despite a phenetic tree is not the 
most important output of a DNA barcoding analysis, it 
can be used to show clearly the pattern of MOTUs compo- 
sition. For this reason, we provided a NJ tree for dataset B 
with MOTUs delimited by squared brackets (see addi- 
tional file 3: 'NJ tree 1 ). 

DNA barcoding performance: comparison of different 
markers and different data handling 

Type C datasets were built with sequences deriving from 
the same specimens in order to allow the comparison of 
the performances of different markers and different data 
managements. The performances were evaluated calculat- 
ing the minimum cumulative error (MCE) rate relative to 
both optimum threshold and Standard threshold values. 
The lesser the MCE, the better is the performance. 
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Figure 2 

K2P distance graph of cox/ filarioid nematodes. Frequency distribution of intraspecific and interspecific genetic diver- 
gences in morphologically identified filarioid nematodes. Graph shows 877 intraspecific and 21775 interspecific comparisons 
across 46 filarioid species. Distances were generated after alignment with MUSCLE, and calculated with MEGA (pairwise dele- 
tion), using Kimuras two parameter substitution model. 


The multiple alignment of coxi gene sequences (dataset 
type C) presents no indels. coxi mean nucleotide distance 
within species is 0.5% (Standard error: 0.5%; range: 0 - 
2.0%); coxi mean nucleotide distance between species is 
15.5% (Standard error: 3.7%; range: 0.2 - 27.8%); coxi 
overall mean diversity is 15.0% (Standard error: 1.0%). 

As expected, the multiple alignment of 12S rDNA gene 
sequences (dataset type C) shows several indels, which 
were in most cases concentrated in the variable regions. 
12S rDNA mean nucleotide distance within species is 
2.2% (Standard error: 1.7%; range: 0 - 6.0%); 12S rDNA 
mean nucleotide distance between species is 17.4% 
(Standard error: 4.2%; range: 0.2 - 34.5%); 12S rDNA 
overall mean diversity is 17.0% (Standard error: 1.1%). 

Based on MCE (relative to OT) rate comparison, the two 
DNA barcodes used show different performances (mean 
MCE for coxi is 0.3% and mean MCE for 12S rDNA is 0.7%). 

Using the marker coxi, the eight different combinations of 
data handling show the same value of MCE indicating 
that the performance of DNA barcoding with the marker 
coxi is not susceptible to the tested data handling. OT and 


ST assume very similar values, and the rates of MCE rela- 
tive to the different thresholds are the same (see Table 2). 

Differently, using 12S rDNA, the eight different combina- 
tions of data handling show rates of MCE remarkably dif- 
ferent indicating that the performance of DNA barcoding 
with the marker 12S rDNA is very susceptible to different 
data handling. With this marker, OT and ST assume very 
different values, and the performance of DNA barcoding 
with the two thresholds is extremely dissimilar (mean 
MCE relative to OT is 0.7%; mean MCE relative to ST is 
50.5%). 

Interestingly, the two markers show very different man- 
ageability: coxi has revealed to be less susceptible than 12S 
rDNA to changes in alignment algorithm, Software used 
for distance estimation, and gap treatment. The lower 
manageability observed for the marker 12S rDNA is cer- 
tainly caused by the presence of numerous indels. 

Also for datasets C, the errors of DNA barcoding per- 
formed with the better data handling are all attributable 
to false negatives and are relative to the couples of species: 
O. volvulus and O. ochengi ; C. bulboidea and C. longa. 
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Discussion 

On type A dataset a really good discrimination level is 
achievable, with 44 species out of 46 identifiable via DNA 
barcoding. For two couples of species, the interspecific 
divergence is less than the optimal threshold and hence 
two morpho-species of filarioid nematodes are not 
resolved by DNA barcoding approach. These two species 
belong to Onchocerca and Cercopithifilaria genus. Despite 
O. volvulus and O. ochengi are easily identified based on 
morphology and host specificity, their nucleotide diver- 
gence is quite low (mean interspecific divergence 1.9%). 
If O. volvulus infects human patients only in Africa (origi- 
nally) and South America (following the transatlantic 
slave trade) and O. ochengi infects only cattles, the two 
species could derive from a recent speciation event [35]. 
This event could decrease the resolution power of DNA 
barcoding. 

Another putative recent speciation has been proposed for 
two species of Cercopithifilaria genus (C. longa and C. bul- 
boidea), showing a mean interspecific divergence of 0.2%. 
These parasites are restricted to two Japanese mammals 
(Naemorhaedus crispus and Cervus nippon), and a recent 
speciation event has also been hypothesized using both 
molecular and morphological data [27,36]. It should be 
noted that these evolutionary dynamies are often difficult 
to identify as reported in [37]. 

Dataset B (that encompasses all the coxi sequences of filar- 
ioid nematodes available in GenBank) has been used to 
perform DNA barcoding and DNA taxonomy with a tree- 
based method. Coherently with the results obtained with 
dataset A, this phenetic approach shows a ciear separation 
of MOTUs representing separated groups of morpho-spe- 
cies with the exception of O. volvulus-O. ochengi and C. bul- 
boidea-C. longa. Anyway, closely related species could be 
characterized by a certain level of interspecific hybridiza- 
tion, because the reproductive isolation could not be total 
since the very beginning of the natural history of a species. 
These effects are particularly evident in mitochondrial 
gene trees, and represent a serious problem for DNA bar- 
coding (at least in most metazoans, for which mitochon- 
drial markers are widely used). Problems ofthis nature are 
likely to have occurred in the O. volvulus-O. ochengi and C. 
bulboidea-C. longa cases where traditional taxonomy iden- 
tified good species [27,38]. As a consequence, the usage of 
a tree-based method alone for species Identification could 
be dangerous and deceptive. Moreover, in a gene tree, a 
'true' species may be wrongly represented by a para- 
phyletic group of alleles/haplotypes, due to introgression 
or incomplete lineage sorting (see [ 15]). In such cases, the 
gene tree could appear misleading or uninformative about 
the species Identification because of retention, and conse- 
quent randorn sorting, of ancestral polymorphisms. 


It is important to underline that GenBank entries are not 
absolutely free from Identification errors. The results of 
DNA barcoding analyses performed on coxi sequences 
obtained from GenBank (dataset B) do not show such 
type of problem. Ilowever, an example of error is repre- 
sented by the entry [ GenBank: AY462911] identified as 
Litomosoides carinii. This species parasites sciurids in Brazil 
[39] and was described by Travassos in 1916. The conge- 
neric species Litomosoides sigmodontis was described by 
Chandler in 1931, parasites the murid Sigmodon hispidus, 
and is spread worldwide in the laboratories as model spe- 
cies for the studies on filarioses. For some reasons there is 
the tendency to confound these two clearly distinet spe- 
cies, and it is relatively common to observe the erroneous 
name L. carinii used instead of L. sigmodontis for labora- 
tory strains of these filariods. In this context, it should be 
noted that basically all the results on L. sigmodontis pub- 
lished till now are relative to these laboratory strains 
established since 1970s. Here we present a molecular 
Identification of L. sigmodontis directly collected from wild 
hosts. Laboratory strains and wild specimens show no 
molecular differences. 

The five unidentified MOTUs present in dataset B encom- 
pass parasites of three avian hosts, a taxonomic group 
where biodiversity and distribution of filarial nematodes 
are underestimated. As described above, these are cases 
where molecular analysis can help to discover new species 
(DNA taxonomy). 

It must be underlined that DNA taxonomy performed 
with simple molecular data can only suggest the presence 
of potential new species, whose real existence must be cor- 
roborated by integrated approaches [40]. 

Type C datasets reveal that two different markers have 
similar discrimination power, but if coxi shows high man- 
ageability in data handling, the marker 12S rDNA is more 
susceptible to the data handling (especially in gap treat- 
ment). Processing 12S rDNA type C dataset with MUSCLE 
and MEGA (pairwise deletion), DNA barcoding performs 
6.3 time better than using MUSCLE and MEGA (complete 
deletion). In addition, processing 12S rDNA type C data- 
set with MUSCLE it is possible to obtain 0.3% ofMCE (see 
Table 2), whilst using ClustalX, it is possible to obtain 
0.4% of MCE (see Table 2). This is a quite relevant obser- 
vation: the generation of a reliable alignment is a major 
impediment limiting the use of 12S rDNA gene sequences 
for barcoding purposes. For this reason, Chu et al. [41] 
have proposed to use ribosomal DNA sequences for DNA 
barcoding without performing an alignment, showing 
congruence between their approach and a tree reconstruc- 
tion (based on neighbour-joining algorithm). Anyway, 
12S rDNA offers practical benefits: it is much shorter com- 
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pared with coxi, and therefore more likely to be readily 
amplified from chemically damaged (i. e. formalin fixed) 
or badly conserved specimens [9]. 

It is important to underline that the presence of nuclear 
mitochondrial pseudogenes (numts [42]) could intro- 
duce serious ambiguity into DNA barcoding and their 
presence cannot be known a priori [43]. In nematodes, 
numts seems to be rare [42], despite their presence has 
been reported (see for example [44] were a short fragment 
of the mitochondrial 1 6S rDNA of W. bancrofti included 
into the nuclear LDR region is used for the screening of 
this parasite). In our study, the results of BLAST search, 
multiple alignment analyses and the quality of trace files 
for bidirectional processing of our sequences seems to 
exclude any interference caused by numts. 

Our results indicate that the proposal to use the ST (10 
times intraspecific variability) as described in [11] must 
be evaluated case by case. Indeed, in the case of coxi, the 
OT is equivalent to ST (both the thresholds generate the 
same value of MCE), but for 12S rDNA OT performs 
extremely better than ST (mean MCE relative to OT is 
0.7%, mean MCE relative to ST is 50.5%). The extremely 
high values of MCE relative to ST are caused by the mod- 
erately high intraspecific K2P distances of the marker 12S 
rDNA that are enhanced of a 10 times magnitude. The 
data handling has also a relevant effect on the mean 
intraspecific divergence: MUSCLE, TREECON and consid- 
ering gaps are all alternatives that enhance K2P distances. 

The sampling of filarioid nematodes is clearly not exhaus- 
tive and particularly difficult, due to complications associ- 
ated with their collection (i. e. recovery at necropsy in 
most of the cases), that requires highly skilled personnel 
and enduring logistic efforts all over the world. The data- 
sets presented encompass also species for which only one 
sequence is available. This is a circumstance that avoid to 
evaluate the intraspecific variability of the marker, and 
consequently the discrimination power of the method 
decreases. However, we want to remark the importance of 
the datasets here reported: filarioid nematodes represent a 
relevant neglected, vector-borne, tropical diseases. 

Conclusion 

DNA barcoding represents a powerful tool for taxonomy, 
but without the integration of traditional approaches 
could become a simple collection of MOTUs. Recent stud- 
ies showed that different approaches to species recogni- 
tion can generate similar results, encountering the favour 
of scientific community [45] suggesting that an integrated 
approach to species recognition is a possibility [46]. In 
our opinion, the establishment, improvement and main- 
tenance of DNA barcoding as a taxonomic tool will 
require a long-lasting interaction between traditional tax- 


onomy and DNA-based approaches. In this work tradi- 
tional and molecular approaches have been considered as 
an integrated method for achieving the goal of species 
Identification. 

DNA barcoding is a good method for taxonomical Identi- 
fication of filarioid nematodes, and it has shown a high 
coherence with classical taxonomy. The results of the inte- 
grated approach to species Identification clearly show 
where DNA-based and morphological identifications are 
consistent, and where they are not. 

This study suggests that both coxi and 12S rDNA appear to 
be appropriate molecular markers for Identification of 
filarioid nematodes at species level via DNA barcoding. 
More in detail, the results of DNA barcoding has been 
shown to be more consistent under different data han- 
dling when performed with coxi than 12S rDNA. On the 
opposite, 12S rDNA is less manageable, but it is easier to 
amplify than coxi. 

The threshold value proposed by Hebert (10 times mean 
intraspecific divergence) [ 1 ] has revealed to be applicable 
for coxi, and not for 12S rDNA. In the case of coxi, the 
threshold value 4.8% can also be used to separate poten- 
tially new filarioid species. We conclude that nucleotide 
sequences of coxi from filarioids are of high interest for 
species Identification throughout DNA barcoding. 
Despite the databases here reported encompass only few 
specimens of medical interests, they represent a useful 
starting point for rapid Identification of these parasites 
and for applications such as epidemiological surveys and 
populational dynamies. 

Using coxi with a threshold that minimise the error rate, 
all of the filarial nematodes involved in medical or veter- 
inarian context (with the exception of O. volvulus and O. 
ochengi) can be coherently identified as morphological 
(species) and molecular entities (MOTUs). Finally, it 
should be noted that the two approaches for species Iden- 
tification (i. e. morphological and molecular) are not con- 
sistent at 100%. This is not unexpected, both methods are 
susceptible to different bias: sibling species, morphologi- 
cal polymorphisms, introgression and coalescence can, 
for instance, confound identifications. One way to seri- 
ously cope with these difficulties is to follow the cross 
control given by different approaches (for works about 
integrated taxonomy see [8,40]). The possibility to iden- 
tify how and where the different approaches are not 
coherent can be the first step for developing of a tme inte- 
grated approach to taxonomy. 

Competing interests 

The authors declare that they have no competing interests. 


Page 10 of 12 

(page number not for citation purposes) 



Frontiers in Zoology 2009, 6:1 


http://www.frontiersinzoology.eom/content/6/1/1 


Authors' contributions 

EF is a post-doctoral researcher interested in bioinformat- 
ics and developed the analyses. MB and AG performed 
PCRs and sequencing. OB, SU, RG and HF collected and 
identified most of the samples used in the analyses using 
traditional taxonomy. CM and CB highly contributed to 
the discussion of the results. CB was also essential in the 
initial contact between the researchers that originated this 
international research network. MC played the role of 
inspirer and coordinator of the research. 

Additional material 


Acknowledgements 

The authors are indebted with the many colleagues involved in the collec- 
tion of biological samples all around the world. We would like to thank Dr. 
Angelica Crottini for critical discussion of the results and two anonymous 
referees for their helpful comments. The molecular work was partly sup- 
ported by the Italian Ministero dell' Istruzione, delTUniversita e della 
Ricerca. 

References 

1 . Hebert PDN, Cywinska A, Ball SL, deWaard JR: Biological identifi- 
cations through DNA barcodes. Proc Roy Soc Lond B Biol Sci 2003, 
270:313-321. 

2. Hebert PDN, Penton EH, Burns JM, Janzen DH, Hallwachs W: Ten 

species in one: DNA barcoding reveals cryptic species in the 
neotropical skipper butterfly Astraptes fulgerator. Proc Nati 
Acad Sci Unit States Am 2004, 1 0 1 : 1 48 1 2- 1 48 1 7. 


3. Blaxter M: Counting angels with DNA. Nature 2003, 
421:122-124. 

4. Floyd R, Abebe E, Papert A, Blaxter M: Molecular barcodes for soil 
nematode Identification. Mol Ecol 2002, 11:839-850. 

5. Markmann M, Tautz D: Reverse taxonomy: an approach 
towards determining the diversity of meiobenthic organisms 
based on ribosomal RNA signature sequences. Phil Trans Biol 
Sci 2005, 360:1917-1924. 

6. Monaghan MT, Balke M, Gregory TR, Vogler AP: DNA-based spe- 
cies delineation in tropical beetles using mitochondrial and 
nuclear markers. Phil Trans Biol Sci 2005, 360: 1 925- 1933. 

7. Vences M, Thomas M, Bonett RM, Vieites DR: Deciphering 
amphibian diversity through DNA barcoding: chances and 
challenges. Philos Trans Biol Sci 2005, 360: 1 859- 1 868. 

8. Lefebure T, Douady CJ, Gouy M, Gibert J: Relationship between 
morphological taxonomy and molecular divergence within 
Crustacea: proposal of a molecular threshold to help species 
delimitation. Mol Phylogenet Evol 2006, 40:435-447. 

9. Casiraghi M, Bain O, Guerrero R, Martin C, Pocacqua V, Gardner SL, 
Franceschi A, Bandi C: Mapping the presence of Wolbachia pip- 
ientis on the phylogeny of filarial nematodes: evidence for 
symbiont loss during evolution. Int J Parasito I 2004, 34:191 -203. 

10. Rach J, Desalle R, Sarkar IN, Schierwater B, Hadrys H: Character- 
based DNA barcoding allows discrimination of genera, spe- 
cies and populations in Odonata. Proc Roy Soc Lond B Biol Sci 
2008, 275:237-247. 

I I. Hebert PDN, Stoeckle MY, Zemlak TS, Francis CM: Identification 
of birds through DNA barcodes. PLoS Biology 2004, 2:e3 12. 

12. Moritz C, Cicero C: DNA barcoding: Promise and pitfalls. PLoS 
Biology 2004, 2:e354. 

13. Meyer CP, Paulay G: DNA barcoding: error rates based on 
comprehensive sampling. PLoS Biology 2005, 3:e422. 

14. Wiemers M, Fiedler K: Does the DNA barcoding gap exist? - a 
case study in blue butterflies (Lepidoptera: Lycaenidae). 
Front Zool 2007, 4:8. 

15. Funk DJ, Omland KE: Species-level paraphyly and polyphyly: 
frequency, causes, and consequences, with insights from ani- 
mal mitochondrial DNA. Annu Rev Ecol Evol Systemat 2003, 
34:397-423. 

16. Besansky NJ, Severson DW, Ferdig MT: DNA barcoding of para- 
sites and invertebrate disease vectors: what you don't know 
can hurt you. Trends Parasitol 2003, 1 9:545-546. 

17. Powers TO: Nematode molecular diagnosties: From bands to 
barcodes. Annu Rev Phytopathology 2004, 42:367-385. 

1 8. Blaxter M, Mann J, Chapman T, Thomas F, Whitton C, Floyd R, Abebe 
E: Defining operational taxonomic units using DNA barcode 
data. Phil Trans Roy Soc Lond B 2005, 360: 1935-1 943. 

19. Pion SD, Clarke P, FilipeJA, Kamgno J, Gardon J, Basanez MG, Bouss- 
inesq M: Co-infection with Onchocerca volvulus and Loa I oa 
microfilariae in Central Cameroon: are these two species 
interacting? Parasitology 2006, 132:843-854. 

20. Anderson RC, Bain O: Keys to genera of the order Spirurida. 
Part 3. Diplotriaenoidea, Aproctoidea and Filarioidea. In CIH 
keys to the nematode parasites of vertebrates Volume 3. Farnham Royal: 
Commonwealth Agricultural Bureaux; 1976:59-1 16. 

21. Stein LD, Bao Z, Blasiar D, Blumenthal T, Brent MR, Chen N, Chin- 
walla A, Clarke L, Clee C, Coghlan A, Coulson A, D'Eustachio P, Fitch 
DHA, Fulton LA, Fulton RE, Griffiths-Jones S, Harris TW, Hillier 
LDW, Kamath R, Kuwabara PE, Mardis ER, Marra MA, Miner TL, Minx 
P, Mullikin JC, Plumb RW, Rogers J, Schein JE, Sohrmann M, Spieth J, 
Stajich JE, Wei C, Willey D, Wilson RK, Durbin R, Waterston RH: 
The genome sequence of Caenorhabditis briggsae: a platform 
for comparative genomies. PLOS Biology 2003, 1 : 1 66- 1 92. 

22. Casiraghi M, Anderson TJC, Bandi C, Bazzocchi C, Genchi C: A phy- 
logenetic analysis of filarial nematodes: comparison with the 
phylogeny of Wolbachia endosymbionts. Parasitology 200 1 , 
122:93-103. 

23. Anderson RC: Nematode parasites of Vertebrates - Their development 
and transmission Wallingford: CAB International; 2000. 

24. Chabaud AG: Le genre Dipetalonema Diesing, 1861; esssai de 
classification. Annales de Parasitologie Humaine et Comparee 1952, 
27:250-285. 

25. Bartlett CM: Cercopithifilaria leporinus n. sp. (Nematoda: Filar- 
ioidea) from the snowshoe hare (Lepus americanus Erxleben) 
(Lagomorpha) in Canada. Annales de Parasitologie Humaine et 
Comparee 1983, 58:275-283. 


Additional File 1 

Investigated specimens. List of species including biological data, acces- 
sion numbers and datasets (for data analysis) ofthe nematodes 
included in this study. Where avaUable date, place ofcollection and hosts 
are indicated (n.d.: no data avaUable). * The host species indicated as 
Naemorhedus crispus is synonymous to Capricornis crispus. ** Labo- 
ratory strain in European laboratories since 1970s. *** Collected from a 
patient travelling from Camerun. **** Collected from a patient travel- 
lingfrom India. 

Click here for file 

[http:/ /www. biomedcentral.com/content/supplementary/ 1742- 
9994-6-l-Sl.pdf] 


Additional File 2 

Experimental conditions. Detailed conditions ofDNA extraction, 
primers, PCRs and DNA sequencing. 

Click here for file 

[http:/ /www.biomedcentral.com/content/supplementary/ 1742- 
9994-6-l-S2.pdf] 

Additional File 3 

NJ Tree. Neighbour joining tree based on coxi sequences generated using 
MEGA 4.0 ( Tamura et al, 2007) - gaps treated as 'complete deletion'. 
Click here for file 

[http:/ /www.biomedcentral.com/content/supplementary/ 1742- 
9994-6-l-S3.pdf] 


Page 11 of 12 

(page number not for citation purposes) 



Frontiers in Zoology 2009, 6:1 


http://www.frontiersinzoology.eom/content/6/1/1 


26. Guerrero R, Martin C, Gardner SL, Bain O: New and known spe- 
cies of Litomosoides (Nematoda: Filarioidea): important 
adult and larval characters and taxonomic changes. Comp Par- 
os/to/ 2002, 69:177-195. 

27. Uni S, Suzuki Y, Baba M, Mitani N, Takaoka H, Katsumi A, Bain O: 

Coexistence of five Cercopithifilaria species in the Japanese 
rupricaprine bovid, Capricornis crispus. Parasite 2001, 
8:197-213. 

28. Kimura M: A simple method for estimating evolutionary rates 
of base substitutions through comparative studies of nucle- 
otide sequences. J Mol Evol 1 980, 1 6: 1 I I - 1 20. 

29. Tamura K, DudleyJ, Nei M, Kumar S: MEGA4: Molecuiar Evolu- 
tionary Genetics Analysis (MEGA) Software version 4.0. Mol 
Biol Evol 2007, 24:1596-1599. 

30. Britten RJ, Rowen L, Williams J, Cameron RA: Majority of diver- 
gence between closely related DNA samples is due to indels. 
Proc Nati Acad Sci Unit States Am 2003, 1 00:466 1 -4665. 

31. Edgar RC: MUSCLE: multiple sequence alignment with high 
accuracy and high throughput. Nucleic Acids Res 2004, 
32:1792-1797. 

32. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The 

ClustalX Windows interface: flexible strategies for multiple 
sequence alignment aided by quality analysis tools. Nucleic 
Acids Res 1997, 24:4876-4882. 

33. Hali TA: BioEdit: a user-friendly biological sequence align- 
ment editor and analysis program for Windows 95/98/NT. 

Nucleic Acids Symp 1 999, 4 1 :95-98. 

34. Peer Y Van de, De Wachter R: TREECON for Windows: a Soft- 
ware package for the construction and drawing of evolution- 
ary trees for the Microsoft Windows environment. CABIOS 
1994, 10:569-570. 

35. Morales-Hojas R, Cheke RA, Post RJ: Molecuiar systematies of 
five Onchocerca species (Nematoda: Filarioidea) including 
the human parasite, O. volvulus, suggest sympatric specia- 
tion. J Helminthol 2006, 80:281-290. 

36. Agatsuma T, Iwagami M, Uni S, Takaoka H, Katsumi A, Kimura E, Bain 

O: Molecuiar phylogenetic relationships among seven Japa- 
nese species of Cercopithifilaria. Parasito/ Int 2005, 54: 1 95- 1 99. 

37. Dorris M, Viney ME, Blaxter ML: Molecuiar phylogenetic analysis 
of the genus Strongyloides and related nematodes. Int J Para- 
sito/ 2002, 32:1 507-1 517. 

38. Bain O: Le genre Onchocerca: hypotheses sur son evolution et 
cie dichotomique des especes. Annales de Parasitologie Humaine 
et Comparee 1981, 56:503-526. 

39. Bain O, Petit G, Diagne M: Etude de quelques Litomosoides par- 
asites de rongeurs; consequences taxonomiques. Annales de 
Parasitologie Humaine et Comparee 1 989, 64:268-289. 

40. DeSalle R, Egan MG, Siddall M: The unholy trinity: taxonomy, 
species delimitation and DNA barcoding. Phil Trans Roy Soc 
Lond B 2005, 360:1905-1916. 

4 1 . Chu KH, Li CP, Qi J: Ribosomal RNA as molecuiar barcodes: a 
simple correlation analysis without sequence alignment. Bio- 
informaties 2006, 22 : 1 690- 1701. 

42. Bensasson D, Zhang D, Hartl DL, Hewitt GM: Mitochondrial pseu- 
dogenes: evolution's misplaced witnesses. Trends Ecol Evol 
2001, 16:314-321. 

43. Song H, BuhayJE, Whiting MF, Crandall KA: Many species in one: 
DNA barcoding overestimates the number of species when 
nuclear mitochondrial pseudogenes are coamplified. Proc 
Nati Acad Sci Unit States Am 2008, 1 05: 1 3486-9 1 . 

44. Abbasi I, Hamburger J, Githure J, Ochola JJ, Agure R, Koech DK, 
Ramzy R, Gad A, Williams SA: Deection of Wuchereria bancrofti 
DNA in patients' sputum by the polymerase chain reaction. 
Trans R Soc Trop Med Hyg 1 996, 90:53 1 -32. 

45. Miller SE: DNA barcoding and the renaissance of taxonomy. 
Proc Nati Acad Sci Unit States Am 2007, 1 04:4775-4776. 

46. Dettman JR, Jacobson DJ, Turner E, Pringle A, Taylor JW: Reproduc- 
tive isolation and phylogenetic divergence in Neurospora: 
comparing methods of species recognition in a model 
eukaryote. Evolution 2003, 57:2721-2741. 


Publishwith Bio Med Central andevery 
scientist can read your work free of charge 

"BioMed Central will be themost significant development for 
disseminating the results of biomedical research in our lifetime . " 
Sir Paul Nurse, Cancer Research UK 


Your research papers will be: 

• available free of charge to the entire biomedical community 

• peer reviewed and published immediately upon acceptance 

• cited in PubMed and archived on PubMed Central 

• yours — you keep the Copyright 


Submit your manuscript here: 
http://www.biomedcentral.com/info/publishing_adv.asp 


O 


BioMedcentral 


Page 12 of 12 

(page number not for citation purposes) 


